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Abstract 

A superparamagnetic nanoparticle (SPN) is a nanometre-sized piece 
of a material that would, in bulk, be a permanent magnet. In the 
SPN the individual atomic spins are aligned via Pauli effects into a 
single giant moment that has easy orientations set by shape or mag- 
netocrystalline anisotropy. Above a size-dependent blocking tempera- 
ture Tb{V, 

Tobs) ) thermal fluctuations destroy the average moment by 
flipping the giant spin between easy orientations at a rate that is rapid 
on the scale of the observation time t^^^. We show that, depite the 
vanising of the average moment, two SPNs experience a net attractive 
force of magnetic origin, analogous to the van der Waals force between 
molecules that lack a permanent electric dipole. This could be rele- 
vant for ferrofluids, for the clumping of SPNs used for drug delivery, 
and for ultra-dense magnetic recording media. 

1 Introduction 

In many areas of physics, forces are effectively suppressed in the interac- 
tion between separated fragments of matter, because of the neutrality of 
each fragment with respect to the appropriate charge quantity. Nevertheless 
"residual" forces still occur between these fragments, typically with a decay 
(as a function of the spatial separation D between the fragments) that is 
different from that of the "bare" interaction. 
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For example, ordinary matter consisting of atoms and molecules is typi- 
cally neutral with respect to electrical charge, but two well-separated charge- 
neutral fragments always experience at least the van der Waals or dispersion 
interaction. This is a residual force that arises because the zero-point mo- 
tions of the electrons on the two fragments are correlated via the Coulomb 
interaction, leading to a non-zero time-averaged force of Coulombic origin, 
despite overall charge neutrality of each fragment. For neutral molecules 
distant D, this leads to an interaction energy varying as -D-^ This is to 
be compared with the bare Coulomb interaction proportional to QiQ2D~^ 
that acts between between fragments with nonzero electric charges Qi, Q2- 
{D^^ is replaced by D^^ when D is large enough that retardation of the 
electromagnetic interaction needs to be considered [I1I2])- 

Similarly the nuclear force between two nucleons has sometines been re- 
garded as a residual color interaction between color-neutral objects. 

Here we propose a similar residual force, of magnetic dipolar origin, act- 
ing between two "superparamagnetic nanoparticles (SPNs)". By this we 
mean that each nanometre-sized particle is composed of a material that is 
ferromagnetic in its bulk state [31 Typically at the temperatures of in- 
terest, the elementary electron spins inside an individual nanoparticle re- 
main locked together by the microscopic exchange interaction, yielding ef- 
fectively a single giant spin with a magnetic dipole moment d^. If the di- 
rections of the giant moments remain steady over time, two such nanoparti- 
cles experience a conventional magnetic dipole-dipole energy proportional to 
dly^do'^ f{6i,(j)i]62,(l)2)/R^- Here R is the spatial separation of the nanopar- 
ticles, and / is a dimensionless function of the angles between each fixed 
moment and the vector R joining the spatial locations of the nanoparticles. 
However each particle has one or more "easy axes" in directions determined 
by magnetocrystalline or, more typically, shape anisotropy. The latter effect 
arises in the strong angular dependence of the magnetostatic self-energy of 
a non-spherical magnetised particle. We will consider the simplest case, in 
which the particle is sufficiently elongated that it has a single easy axis, i.e. 
dominant uniaxial shape anisotropy. Then the energy of a single nanoparticle 
is lowest when its giant spin (dipole moment, d) lies parallel or antiparallel 
to this easy axis. Because the energy barrier Eo for rotation of d between 
easy orientations (not mechanical rotation of the particle) derives from the 
magnetic self-energy of the nanoparticle, it decreases with decreasing volume 
of the nanoparticle. For very small particles, therefore, the projection of d on 
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a measurement axis averages to zero over time, because of repeated thermal 
flipping of the giant spin [3] , caused by thermal agitation from the heat bath 
(e.g. a fluid or solid matrix) that surrounds the nanoparticle. Thus on time 
average the nanoparticle is "neutral" i.e. it has a zero magnetic moment. 

When the thermal agitation of the giant spin is insufficient to flip it be- 
tween easy orientations within the observation time, Tabs, the nanoparticle 
is "blocked", i.e. apparently frozen as to its magnetism. This occurs below 
the blocking temperature of this nanoparticle, Tb, which depends on Eo and 
therefore on the volume V of the nanoparticle. If the relaxation time of d 
over the barrier £'0 is r, then Tb is defined by t(T{,) = Tabs- Blocking is thus 
a purely dynamic phenomenon: extending the observation time, or lowing 
the frequency, lowers Tb and vice- versa [3J. For the present case of SPNs 
suspended in a fluid, the observation time Tabs will be a relevant time for 
mechanical motion of the SPN through the fluid. - e.g. a rotational or trans- 
lational diffusion time. Note that the direct dipolar magnetic interaction 
between SPNs could in principle lead them to clump. However when T > Tb 
the motion of the SPNs through the fluid will not " see" the bare dipolar mag- 
netic interaction between the SPNs, as it has been averaged away between 
attractive and repulsive values during the thermal flipping of the spins. It 
could lead to additional Brownian type of damping and difusion of course, 
but we show here that there is also a net attractive force between SPNs even 
above the blocking temperature. 

The destruction of the permanent magnetic moments by thermal flucta- 
tions is highly undesirable in the case of a magnetic data recording medium, 
where very fine magnetic particles in the nanometer size range will be needed 
in order to pack the magnetically stored data as densely as possible for the 
next generation of devices. The thermal destruction of the permanent mo- 
ments means that data cannot be stored over long times. 

On the other hand, as will be discussed below, the same thermal flipping 
occur ing for T > T;, is beneficial in the case of nanoparticles deliberately 
suspended in human blood as carriers for drug or thermal therapies, since 
now the clumping of the nanoparticles from magnetic dipole interaction is 
suppressed because each particle has effectively a zero magnetic moment. 
The strong clumping that would occur for fully ferromagnetic particles from 
their dipole-dipole interactions could be clinically dangerous, poten- 
tially causing blockage of blood vessels, difficulty of elimination etc. We will 
show below, however, that despite the vanishing of the average individual 
moments, there is a residual attractive interaction between two superparam- 
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agnetic nanoparticles separated by distance R, that falls off as (const) / . It 
is the magnetic analog of the van der Waals or dispersion force that arises via 
the Coulomb interaction between fluctuating electric dipoles on two electri- 
cally neutral molecules [5] lacking permanent dipole moments. This residual 
force could also lead to clumping of the nanoparticles, and so its analysis 
could be signifiant in modern magnetic-particle therapies 



2 Simple preliminary model 

The model described here is based on an argument frequently used to explain 
the attractive van der Waals energy proportional to —R~^ that arises between 
temporary electric dipoles occurring on a pair of electronically neutral atoms 
separated by distance R (see e.g. [5]). It is not rigorous derivation, but may 
help to elucidate the more careful and general mathematical treatment to be 
provided in later Sections. Consider two superparamagnetic nanoparticles 
SPNl and SPN2 as defined above. While averaging to zero over time as 
described above, the magnetic moment on SPNl can exhibit a short-lived 
thermal (or quantal) fluctuation so that its value S^\t) is nonzero at some 
particular time t. For simplicity we will assume that only magnetizations of 
SPNl and SPN2 along one axis (say z) are possible so that S^\t) = S^\t)z, 
and we will consider the case that the spatial separation R between SPNl 
and SPN2 is parallel to x. Then the spontaneous moment di{t) produces a 
dipolar magnetic induction (B-field) 

p)(^t) = -^R-'S'\t)z 
Atx 

at the position of SPN2. Responding to this field, SPN2 produces its own 
magnetic moment 

3^\t) = x^%^^\t) = -x^^)l^R-^S^\t)z (1) 

4:71 

where x*^^"* is the dynamic magnetic susceptibility of SPN2, assumed for now 
to represent an instantaneous response to the field. (Note that here x rep- 
resents the response of the total magnetic moment of the SPN to a small 
applied magnetic induction b. By contrast, the symbol x is normally used 
for the response of the magnetic moment per unit volume to a small applied 
magnetic field h. Thus for a single SPN of volume V, 

d = xb; ^ = xh; X= ^ (2) 
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The dipole ([T]) in turn produces a dipolar magnetic induction back at the 
position of SPNl: 

The interaction energy of this back-field with the original moment d^^\t) is 

E = -6«(t).c?i)(t) 
and this energy has a time or thermal ensemble average 

which is non-zero because (^{S^\t))'^'^ ^ even though {S^\t)) = 0. 

This negative energy produces, upon differentiation with respect to R, a 
net time-averaged attractive force between SPNl and SPN2 that falls off as 

The above simplified theory produces the basic physics and the force, 
but it glosses over a number of issues, such as the role of entropic effects at fi- 
nite temperature, the tensor nature of the magnetic dipole-dipole interaction, 
the quantal aspects of the problem, and the retardation of the electromag- 
netic field. Also, the response X2 has been assumed to be instantaneous, 
whereas there can be a strong and important frequency dependence (time- 
delayed aspect) to the linear response of a SPN. All of these considerations 
are treated in detail in the theory given the next Section. 

3 Detailed theory 

The magnetic dipolar energy (hamiltonian) between two particles with mag- 
netic dipoles dJ-^^ and separated in space by a nonzero vector R = RR, 
is of form 

^(12) ^ _^-3 ^ df%{R)df (3) 

ij=l 

where 

rjn _ f^O 3-Ri-Rj — SijR 
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We assume that we are above the blocking temperature, T > Tj,, i.e. that 
the temperature is high enough (compared with the anisotropy energy bar- 
rier), that each isolated giant magnetic dipole has zero thermal expectation 
taken over the time-scale of interest 

< Jti) >o = = < >o . 

The theory to be developed here is meaningful provided that the thermal 
fluctuations of the moment occur on a time-scale r that is short compared to 
the time Tobs = T^ech ^or the nanoparticle to change its spatial position (or 
physical angular orientation) appreciably, within its fluid medium. Under 
these conditions we will derive a residual attractive force between the two 
superparamagnetic nanoparticles, that could for example be used to study 
residual clumping effects in fluid suspension at temperatures above the block- 
ing temperature . 

The quantum-thermal expectation, denoted < >, of the interaction en- 
ergy between the giant spins is 

^(12) ^ ^ ^(12) ^ ^ _^-3 y^^.(^) < ^(1)^(2) > 

ij=l 

However at flnite temperature it is not this energy but the corresponding 
thermal Helmholz free energy 

A = < > -TS 

that must be considered, where S is the entropy. We achieve an expression 
for A via a Feynman-theorem argument in Appendix A for a classical treat- 
ment of the fluctuations, and in Appendix B for the fully quantal case. In 
either case the result is 

^(A = 1, T, R) - A(X = 0, T, R) = /' (4) 

Jo A 

Here the subscript A means that the quantity is evaluated in the thermal 
ensemble with modified interaction 

H^^^^ = A//(i2) = -XR-^ df%{R)df\ < A < 1 (5) 

E^' Hf^ >, . (6) 
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Since the coupling will be zero (equivalent to A = 0) at infinite separation 
i? — i> oo, we can write Eq (jl]) as an expression for just the free energy of 
interaction between the two nanoparticles: 



A{T,R)-A{T,R^oo) = f Ef)-^ 

JO A 

= -R-^ T,,{R) J < dfUf >x dX (7) 



The problem now reduces to the calculation of the equal-time cross-correlation 
function < df'^d^p >x between the moments in a thermal ensemble with A- 
reduced interaction. 

The equal-time correlation function < d^l^df^ >x can be recovered from 
the time Fourier transform 



gx{uj) = Gx{t)exp{-iujt)dt = / < d]^ {0)d) ' (t) >x exp{-iujt)dt 

of the time-displaced correlation function Gx, 

Gxit)^<d^\0)df\t)>x . 

We can use the finite-temperature fluctuation-dissipation theorem (see 
e.g. [6]) to relate the fluctuation quantity < d^^^df^ > to the dipole-dipole 

(12) 

response function xlj of the combined interacting system, defined in Eq 
([TOD below: 

< d^^df >x,. + < dfdf^ >x,.= ^_ Jp^^^^^ /m {xSj'V + zO) + xj- W + ^0)} 
where /? = -j^. Then 



/-I \ (rf\ /o\ /-I \ X 

< d\ 'd) + d) 'd\ >x,equaitime= — / 5'a(^) exp(iu;t) 

ZTT J-oo 



-OO 7r(l — eyi\i{[3huj)) 



^ Im (c. + .0) + xJS i^ + ^0)}du; 



This can also be expressed as a Matsubara sum by closing upwards in the 
complex u plane, using Cauchy's theorem to obtain a sum of residues at the 
poles Un = iUn = i2n7r / {jSk) , but we will not make explicit use of this here. 
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The interaction energy E and Helmholtz free energy A then become 

3 

7r(l — exp{(3huj)) 



>l(fl)-A(oo) = /^'Bfl ^ (9) 



We assume we know the dipole responses Xij^i^)^ Xij\{^) of each iso- 
lated giant dipole SPNl, SPN2 to an external B field b such that 

d^\Lu) expi-tujt) = ^ x^\uj)b^'^ expi-tujt). 

3 

These individual responses must express the known superparamagnetic prop- 
erties of individual systems. In general it should also describe any Brownian 
tumbling aspects of the response, in the case that the time scale of these 
tumbling motions overlaps that of the magnetic reponse behaviour of each 
SPN. For now we assume that the tumbling is slow so that only the magnetic 
response of a SPN oriented in a fixed spatial orientation is required to appear 
in x.. (The interaction energy may of course depend on the details of this 

— (12) 

orientation, which will be manifested in the particular values of xlj iii the 
chosen cartesian frame.) In Appendix C we discuss a simple model for the x 
of a single isolated SPN. However, to calculate the interaction of two SPNs, 
Eq ([8]) requres knowledge of the cross-response function (cross-susceptibility) 
X^^^ for the interacting pair of SPNs. This is defined as the linear response 
of SPNl's moment to an alternating B field that acts upon SPN2 only: 

^•(^) = x!S(^)^f ' (10) 

where the subscripts i,j label caretesian components of the vectors. To 
calculate x*^^^"* "W'e now consider the slightly more general situation where 
independently-specified small external B fields exp(— icut), exp^—itot) 
are applied to the individual dipoles, in the presence of the dipolar coupling 
between the two systems. 

In time-dependent mean- field theory (RPA), the equations of motion of 
the coupled systems are (all at arbitrary frequency a; and with Einstein sum- 
mation convention for repeated indices): 

= + Ai?-%,rff ) (11) 

df = xg ) + Ai?-^T.,<)) (12) 
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These equations describe the evolution of each giant spin in an effective B 
field containing a time-dependent contribution due to the polarizaton of the 
other giant spin. Using ([12]) to eliminate in(fTTl). we get 

Then for 6^ = (i.e. an external oscillating B field applied only to moment 
SPN2) we have 

^(1) = y (12) 5(2) 

where 

(13) 

c(12) _ ^W\p-3rp ^(2) 
( 1 - \'R-'x^X^xfeTel -\'R-'x^XpxfeTe2 R'^Xl^^xfeTe, \ 

e, = -X'R-'xSX^xfeTs, 1 - \'R-'xSXpX%Ts2 -X'R-'xSX^xfeTe, 

V ->^'R-'xiXpxfe%^ -\'R-'X^X,xTeTe2 1 " \' R-'xi^.X^^TsS J 

This becomes simpler if we have strictly uniaxial responses of the individual 
spins along (say) the x axis, i.e. 

and similarly for x'-^)^ Then we can ignore the 2 and 3 components ^2^'', ^3^^ 
and only need solve a scalar equation, giving 

v(-) = >^x'''R-'TuX^'' (u) 
1 - A2i?-6^(i)Tnx(2)Tn ^ ' 

and then from 

) = -\^R-^T^, r -A^ -Im —J^l^L—^du 

From ([9]), the corresponding free energy of the residual interaction is 
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A{R)-Aioo)= 

Jo A 

The corresponding force between SPNl and SPN2 is 

Eq f|T5|) is valid for the uniaxial case but is readily generalized: there is in gen- 
eral a sum of logarithms of the eigenvalues of the matrix 1 — R~^Tx^^^Tx^'^\ 
Note that both x^^^ and x*^^) in f[T^ \TE\} are frequency-dependent. If the 
denominator of (IT^ vanishes for some frequency Uoj then we have a finite 
oscillation of the magnetic moments for zero driving field - i.e. a free magnon 
collective oscillation mode of the coupled giant spins. Indeed the free energy 
f lT^ can be related to a sum of the thermal free energies of these magnons. 
Actually for the present model, namely x{^) = Xo/(l ~ "^^t) (see Appendix 
C) these frequencies will have a large imaginary part (damping), so there 
are really no magnons in the absence of an applied DC magnetic field. The 
exception is the case ~ 0, where the damping vanishes. If one of the 
magnon frequencies vanishes, ujoj = 0, then we have an instability and the 
system will try to "feeze in" the magnon. This means that the denominator 
in (HM vanishes for zero frequency, which, as the coupling is increased, will 
happen first for A = 1, i.e. 

l-i?-^X^i)(0)x(^)(OKi = 0. (17) 

In time-dependent mean-field theories such as this, this behaviour is usually 
taken to indicate a transition to a broken-symmetry state - in this case the 
moments presumably freeze into a permanent ordering in the antiparallel 
configuration. 



4 Energy in second order (weak coupling) 

Note that if we only want the energy to second order in the interaction then 
from ([8]) we only need x*^^^"* to fii'st order in Tjj, so we can take = in 



10 



CS]), giving 



xlS-AxSi?-%.xl?. (18) 



1 \2dA _ 1 



SO that (jH]) becomes, since /q A^ — 



2' 



residual 



A(R) - Moo) = 1^ Er-'^-^ (19) 

1 ^ 

--R-^ ^ T,j{R)h,^,, where (20) 



— 7r(exp(/3aa;)-l) 



Xll^iuj + tO)T^,{R)x]:Jiu; + tO)\ du. (21) 
The dependence is apparent. In the case of uniaxial response 



A'"'"^''^^ = -^R-^Tii{R)Tn{R)I, where (22) 



h 



n{exp(phuj) — 1] 



Im 



X^^^ {to + iO)x^^^ {uj + iO)] dw (23) 



From Appendix C, a simple model for a superparamagnetic susceptibility 
is x{^) = Xo/(l —iujr), and the frequency integral I in ( 122|) can be estimated 
analytically in two limits depending on the thermal flipping time r of the 
giant spins (see Eqs (136!) and (!35|) of Appendix C). This gives a residual free 
energy 

^residual ^ d-6t-i2 -2 J ^B^; ^bT" >> hr ^ , . 

A =--R TnXoji^^-i^ kBT«hr-' ^^4) 

and a residual van-der-Waals-like force 



Q /I residual ( U 'T 

riresidual o D-7m2 -2 J '^B-' ) 



7t^2 -2 / ^B^, » hr ^ 

kuT « hr-^ 



(25) 



5 Orders of magnitude 

5.1 SPNs below the blocking temperature 

First consider the energy and force of interaction between two SPNs below 
their blocking temperature so that each has a permanant magnetic moment 
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of magnitude = rijiB- At separation R the direct dipole-dipole energy is 
dependent on orientation but is of order 



j^direct 



^dlR-' = (10-^)(9 X 10-^^)^(10^)=^ i?/10 



An 



^-9 



-3 



-9 N 3 



8.1 X 10~^V ( ^"^^"^ 1 Joule (26) 



For example if n = 1000 and = 1 nm, E'^^''^^* ^ lO'^o J. At T = 300K the 
thermal energy is kBTroom = 4 x lO'^V , so ^ 2fcijT. Thus if the two 

SPNs are not thermally suppressed at T = SOOi^' and are able to approach 
to within a nanometer, they will not be prevented by thermal effects from 
rotating to the antiparallel configuration and binding (clumping). 

The corresponding force _p°'*''^'=* between the SPNs is highly orientation- 
dependent but is of order 



^direct 



^UlR-^ ^ 3 X 8. 1 X 10-2^(10-9)-in2 
47r ^ 




2xlO-iVP^-^ N. (27) 



For n = 1000 and R = 1 nm this gives a force of order 20pN, which is small 
but should be directly detectable via Atomic Force Microscopy (AFM) with 
single SPNs attached to substrate and tip. 

5.1.1 SPNs above the blocking temperature 

Now consider a similar system but with a blocking temperature below room 
temperature so that at SOOi^ there are no permanent moments. Then the 
vdW-like theory derived above gives the free energy of interaction. For 
numerical estimates we assume uniaxial susceptibilites and work in the weak- 
coupling limit. We also assume that the giant spins have a zero-frequency 
susceptibility 



corresponding to a giant moment of n Bohr magnetons. Then gives 

^res^dual^^ '2 X IQ-^ {^fhr-^ J oule, hr'^ » k^T 
-2 X 10-% l^^^y keT Joule, hr'^ « ksT 
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For example, let n = 1000, R = 1 nm , T = 300K, and hr « ksT . Then 
j^residuai ~ _2 X lO-^^^^fc^T = 2kBT. Tliis means that for the present case 
the residual energy predicted by the perturbative theory is about the same as 
the direct energy ( l26l) . which is unphysical and simply means that the weak- 
coupling condition is not met and we need (at least) the full RPA theory here 
(Eq. (fT5|) ). If we are in the limit hr^^ » ksT the residual interaction will 
be even larger. In this case the system of two SPNs, despite the thermal 
averaging of an individual SPN, is most probably near to a trasition to a 
spin-locked configuration. In the RPA theory the onset of this condition 
would correspond to a zero denominator in f|T^ . This would occur for 
R-'yS^TuyS^Tn « 1 , i.e. 

For the present case with n = 1000 and T = SOOi^, the crossover occurs 
at about Riock = 1 nm, which is consistent with the above finding that the 
perturbative calculation of the attraction at this separation was unphysical. 

To give another example, suppose that n = 100, R = lOnm, T = 
300K, and hr'^ « keT .Then A^^^^^^""' ^ -2 x l{)-^^^^{^fkBT = 2 x 
lO^^^/c^^; whereas the direct interaction between permanent moments under 
the same conditions from ([26D is E'^^''^^* = 8.1 x IO^^t^qO^ (^) V(300(1.24 x 
lQ~'^^))kBT = 2 X 10~^A;bT. So for this example, neither the direct nor the 
residual interaction would tend to lock the SPNs into an antiferromagneti- 
cally aligned pair. The mechanical forces on the SPNs due to the spin-spin 
interaction in either the direct or the thermally smeared residual case would 
be negligible in the context of normal Brownian motion. 



6 Clumping considerations in fluid suspen- 
sion of SPNs 

(i) Consider SPNs with n = 1000 and with a bocking temperature satisfying 
Th > SOOi^' . From the numbers shown above, at T = 300i^', if they are able 
to aproach one another within about a nanometer, these particles will form 
pairs or larger clusters ("clumping") that are due to the direct (not thermally 
suppressed) magnetic dipolar interaction, and that are not readily broken by 
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thermal processes. Furthermore under these same conditions the demagneti- 
zation field inside a single SPN might be significant, so that the SPN would 
no longer contain a single domain as assumed so far. At larger separations R 
the binding energy falls off as -R"'^, and so the direct magnetic energy, as R 
is increased, will soon be less than the thermal energy ksT. The interaction 
at these larger separations will not immediately cause binding, but may well 
determine the kinetics of closer approach between nanoparticles, resulting 
ultimately in clumping when shorter separations are attained. This process 
is complicated by the strong orientational dependence of the direct interac- 
tion ([3]). SPNs will tend to rotate mechanically within in the fiuid, in order 
to minimize the free energy in the " antiferromagnetic" relative orientation, 
after which their mutual force is attractive. Thus the kinetics of clumping 
will be far from straightforward. 

If clumping is undesirable, the n'^R~^ dependence of of the direct SPN- 
SPN magnetic binding energy suggests that smaller SPNs (e.g. n = 100) 
will be desirable because they are less susceptible to clumping, i.e. they 
can approach to smaller distances (e.g. R = v^lO~^nm = 0.2 nm) before 
clumping occurs. In fact, at such small separations R, the point dipole 
approximation used here may break down, softening the interaction and pos- 
sibly leading to the conclusion that the binding energy even at contact is less 
than the thermal energy. This would imply minimal clumping. 

(ii) Consider SPNs in suspension at T = 300K, with n = 1000, but now 
with Tfy < 300K. Here, despite the thermal suppression of the net individual 
moments, there is a uniformly attractive residual magnetic SPN-SPN free 
energy _ This varies as n^R~^ within the perturbative approxima- 

tion (see fl29l) ). and so becomes much weaker than the direct interaction at 
large separations R . However at shorter separations, the stronger n and 
R dependence of the perturbative residual energy expression fl2^ suggests 
that E^'^^^'^ could exceed E'^^'^'^^^. This is of course unphysical: the correlations 
between the orientations of giant moments that give rise to E^^^'^'^ cannot be 
greater than perfect correlation, corresponding to the direct interaction in 
the antiferromagnetic configuration of the two giant moments. Thus in gen- 



eral 



E 



resid 



< 



E 



direct, max 



In fact the perturbative approximation breaks 
down in small-i? regime, and the full RPA expression f|T5l) will be needed 
instead of fl29|) . We do not yet have analytic energy and force expressions 
owever it is clear that this approach can yield a residual 



in this regime. I: 



interaction 



resid 



of a strength approaching 



E' 



direct, max 



It seems likely. 
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therefore, that because of the residual interaction, there will not always be 
a discontiunous cessation of clumpimg as the temperature is raised above 
the blocking temeparure Tb. However the direct interaction can be repulsive 
whereas the residual interaction is always attractive, so there is scope for 
some quite rich behaviour. 

7 Prospects for experimental verification of 
the theory 

7.1 Direct measurement of the force between two in- 
dividual SPNs 

In Section 15.11 above, the direct interaction between permanently magne- 
tized SPNs with n = 1000 at separation R = 1 nm was estimated to exceed 
kBTroom, and the force was estimated as 20 pN. A force of this magnitude is 
likely to be observable, with some care, via atomic Force Microscopy. The 
simplest configuration might involve one SPN attached to a non-magnetic 
substrate, and another SPN attached to the AFM tip. One could then 
measure the force as a function of temperature. One might expect a reduc- 
tion in the measured force as T is increased above T^. As discussed above, the 
force could even change from repulsive to attractive, depending on the initial 
orientation of the giant moments prior to heating and subsequent destruction 
of the net moments. The need for a measurably large force puts us out of 
the perturbative regime for the residual interaction, so more straightforward 
but messy theoretical work will be required in order to predict the way in 
which F varies with distance and temperature near {R,T) = {lnm,Th). It 
is not clear whether the force will be large enough for AFM detection in the 
regime of larger separations where the perturbative analysis f l29|) is valid. 

7.2 Indirect measurement via observation of structure 
factors in fluid suspension 

Here we propose (e.g.) small-angle xray diffraction measurements on SPNs 
in suspension in a viscous fluid such as glycerine. The metallic SPNs should 
provide good Xray contrast. The measured structure factor of the array of 
SPNs should reveal evidence of positional corrrelations between the SPNs, 
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which in turn is related to the forces between the SPNs as predicted here. 
Again, one hopes to sec some changes as the temperature is raised through 
the blocking temperature Tb. 

7.3 Magnetic resonance experiments 

Although the present theory did not predict any lightly damped magnons 
(combined oscillations of the magnetic moments) for a pair of adjacent SPNs, 
there might be the possibility of such modes if a strong DC magnetic field 
is applied. Magnetic resonance experiments might then be able to detect 
shifts in the single-SPN resonance frequency due to the proximity of a pair 
of SPNs. Even without the external DC field, an analysis of the linewidth 
of the zero- frequency "resonance" might reveal information about SPN-SPN 
coupling. 

8 Summary and future directions 

We have predicted a residual force between superparamagnetic nanoparticles 
that persists above the blocking temperature. The force is the magnetic ana- 
logue of the electrically-driven van der Waals interaction between electrically 
neutral molecules. Our theory also deals with the dynamic spin response 
of coupled SPNs to small ac external magnetic fields. Our results may 
be experimentally testable, and may have implications for ferro-fluids, for 
nanoparticle-based medical therapies, and for magnetic recording technol- 
ogy. The new force is most likely to be significant for nanoparticles that 
approach one another quite closely, at separations of 0(nm). At these sep- 
arations the point-magnetic-dipole approximation used here will need to be 
replaced by a theory that attributes a finite spatial size and definite physi- 
cal shape to the nanoparticles. A good starting model will be an ellipsoidal 
shape, and fortunately the full electrodynamic theory of Casimir interactions 
is quite well developed for this geometry. A theory along these lines will be 
the next step. 
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9 Appendix A: How to deal with the entropic 
part (classical angle-distribution approach) 

The joint state of two interacting superparamagnets is specified by a clas- 
sical distribution /(^^(f2i,f22) in the two solid angles fii, ^2 defining spatial 
directions where the 2 giant spins point: 

The reduced-strength interaction XE between the superparamagnets is 
given by ([6]). Then from general thermodynamic principles, at a given tem- 
perature T , coupling strength A and separation R, the correct distribution 
f^^ (T, R) is that which minimizes the trial free energy: 

A(A,T,i?) = Mm^(,)^|^(,)|^^A(A,r,i?: 

so that the following functional derivative is zero 

6A 



5p) 



where (30) 



A{X,T,R: = XE-TS 

= XR-^'j2T^J{R) E df\Q^)df\Q2)f^'\n^,Q2) 

ij=l Ui,U2 

+kBT E f'\n^,Q^)\nf^\Q^,Q,) 

Consider an infinitesimal increase in the coupling strength from A to A + AA. 
As a result, f^'^^ changes by an amount A/^^^ and noting that E =< H >\= 
E\/X we have a resulting change in A : 

AA = AAE + E ^ A/(2) = AA:^ + 



where the zero comes from (l30l) . 

Notice that we only have to know the interaction E and not the entropic 
part, to find the change in A. 

Then the change in A in switching on the interaction adiabatically is 



Ax=i-Ax=o= I\a= f 
Jo Jo 



1 dA 
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We have already shown how to calculate Ex by using the fluctuation- dissipation 
theorem and the mean-field (RPA) assumption. Also note that the A = 
value of the free energy is independent of separation R: 

\ni Hz / 

Thus the entire independence of A{X = 1,R,T) is captured by the integral 

10 Appendix B: How to deal with the en- 
tropic part (fully quantal approach) 

Our quantum mechanical basis (NOT the eigenstates) for the combined mag- 
netic state of the two systems together consists of the factorised states 

Ki) = K) \j) ■ 

where the first ket refers to quantum state of SPNl and the second ket to 
SPN2. 

The thermal density matrix operator of a pair of magnetically interacting 
nanoparticles has matrix elements in this basis denoted by 

Pij:kl 

and traces can be taken over this or any other basis with the same result. 

We consider starting from the thermal equilibriumn of two isolated nanopar- 
ticles, and consider the effect on the free energy of turning on the interac- 
tion by replacing the inter-nanoparticle interaction hamiltonian H'^^'^\R) by 
XH^^^\R), and then increasing A from to 1 while holding the inter-particle 
separation it! fixed. 

For coupling strength A the Helmholtz free energy is a trace: 

A{X, T, R, =E-TS = Tr ((^o + XH^''\R)) p^^)) -fc^Tr (p^^) lnp(2)j 
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For fixed Hamiltonian, at thermal equilibrium A is stationary with respect 
to arbitrary variations in density matrix that preserve Trp^'^\ see e.g. []): 

so that 

Then the first-order change in the equilibrium free energy, when the coupling 
is increased from A to AA, is 

= AATr {h^'^\R)p^) + ATr (^Ap^^)^ = AA^^ + 

Then the change in free energy in switching on the interaction between the 
two systems is 

A(A = l,T,i?)-A(A = 0,r,i?)= (32) 

JO A 

The same formula can be derived for the classical case, by considering 
a pair distribution /(f2i,fi2) of angular orientations VL of the two giant mo- 
ments. (See Appendix A). 



11 Appendix C: simple superparamagnetic model 

for Xap(^) 

Assume that the individual giant moment has its easy axis along e. Then 
the response to a field h is only via /i.e, and the response is along e 
i.e. 

Here we assume a widely-used model pi] for the frequency-dependent mag- 
netic susceptibility of an individual SPN in the absence of a d.c. external 
magnetic field: 

X{^) = — ' ^ = Toexp{Eo/{kBT)) 

1 — tUJT 
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where r is the thermal flipping time of the giant moment, assumed to arise 
from an intrinsic attempt time Tq and a thermally-activated Boltzmann suc- 
cess rate in surmounting the anosotropy energy barrier Eq. Then 



Im 



Im 



Im 



-(2) 

Xo 



1 — zcjr^^) 1 — iujT^'^^ 
xS'^(l + z^r«)x?^(l + za;r(2))' 



1 + cj2r(i)2 1 + cu2r(2)2 
^ X?^X?^o^(r(^) + r(2)) 

(1+C^2^W2) (1+^2^(2)2) 

If the two superparamagnetic nanoparticles have the same parameters, we 
get 



Im 



(1 +cj2r2)^ 



(33) 



Then the second order energy fl20|) between two nanoparticles with easy 
axes e^^), e*^2) becomes 



n 



ij=l 



7T{exp{(3fiuj) — 1) 



xlmxl'Jico + iO)R~%,iR)t3iu + tO)di 



(2) 



ij=l 



-x'^^\uj + iO)x^^\uj + iO)du 



'-00 TT{exp{/3?iuj) — 1)' 
For two similar SPNs the frequency integral is 

h 2xluJT 



7r(exp(/5M-l)(l+o^V2)'""' ^^^^ 

If jSh/r << Ithen the Imx^ factor cuts the integral off for \uj\ > t^^ , i.e. 
for (3huj > 1 and we can Taylor-expand the denominator to 1st order giving 



h 



-oo TT/Shu (1 + W2r2)^ 

n r 



TT 



-oo {l+UJ^\T\y 



-doj 



TT 



(35) 
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On the other hand if pU/r » 1 then 



TT J-oo 



duj — 



OJ 



TrduO 



TT 



TT 




-1 



(36) 
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